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Abstract. In a ~ 20 ksec XMM-Newton observation the X-ray transient radio-quiet quasar PG 0844+349 was found in 
a historically high state compared to previous X-ray observations. The quasar showed a featureless spectrum with a strong 
soft excess over the extrapolation of a hard power law. Comptonization models or a broken power law with Fsoft ~ 2.75, 
Fhard ~ 2.25 and a break energy of Ebreak ~ L35 keV represent acceptable descriptions of the spectral continuum. In the 
Comptonization models the temperature of the Comptonizing gas is considerably lower than generally found in (broad line) 
Seyfert galaxies whereas the optical depth is much higher. As a similar behavior has been seen in NLSyl galaxies, it might be 
an indicator of the different physical conditions in these two classes of AGN. During the XMM-Newton observation the flux of 
PG 0844+349 varied achromatically in a smooth, nearly linear fashion, by ~ 25% on time scales of a few thousand seconds, 
which puts some constraints on current models of Comptonizing accretion disk coronae. 
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1. Introduction 

Quasars as a class of Active Galactic Nuclei (AGN) exhibit 
a vast diversity in their observed properties. Some of these, 
like the radio-loud - radio-quiet dichotomy, can be related to 
the physical characteristics of the accreting central black hole 
and to morphological differences of the host galaxies. Others 
can result from specific geometrical conditions of the accretion 
flow and from the particular viewing condition of the observer 
with respect to the quasar (see, for example, Elvis 2000). In 
many cases, however, the causes for the unusual properties of 
a quasar remain obscure and thus the study of objects display¬ 
ing atypical characteristics might yield important clues to the 
physical processes governing the quasar emission. 

One of the objects with several extraordinary properties 
is PG 0844+349, a well studied nearby (z=0.064), bright 
(mv=14.0), radio-quiet quasar from the Palomar Green sample 
(Schmidt & Green 1983). Deep radio observations reveal no 
nuclear source and the faint detection in the field (81,49 GHz = 
0.3 mJy) might be related to an R = 20 mag object south-west 
of the quasar (Condon et al. 1987). PG 0844+349 is optically 
quite spectacular: it appears as a barred spiral with complex 
outer structure and a nearly equally bright companion galaxy 
(Hutchings & Crampton 1990). It belongs to the class of AGN 
with unusually strong Fe II emission (Wang et al. 1996) and ac¬ 
cording to the optical spectrum displayed in Boroson & Green 
(1992) the source shows all features of a Narrow-Line Seyfert 


1 galaxy (NLSyl) - strong Fell emission and weak forbid¬ 
den narrow lines, although its FWHM(H/3) = 2420 km is 
slightly above the often used cut-off line of 2000 km (e.g. 
Osterbrock & Pogge, 1985). However, this cut-off for the clas¬ 
sification of NLSyl s is more artificial than physical, so we still 
might define this source as a NLSyl galaxy especially as it 
shows X-ray properties typical for this class of objects, i.e., a 
steep soft X-ray spectrum and strong variability (for reviews 
see Boiler et al. 2000). It is infrared-loud with a disturbed host 
galaxy (Clements 2000), and it is strongly variable at all wave¬ 
lengths on various time scales. The optical flux seems to have 
changed by ^70% between 1986 (Elvis et al. 1994) and 1993 
(Maoz et al. 1994) and intra-night variations with amplitudes of 
AB — 0.08 mag were detected by Jang et al. (1997). From re¬ 
verberation measurements Kaspi et al. (2000) deduced a mass 
of ^ 2.5 X 10^ M 0 for the quasar. 

Compared to the ROSAT All Sky Survey (Yuan et al. 1998) 
the soft X-ray flux of PG 0844+349 had decreased by a factor 
of ~6 in a pointed ROSAT observation 6 months later, with¬ 
out any noticeable changes in the spectral slope (Rachen et al. 
1996). During the X-ray low state PG 0844+349 had an X- 
ray loudness of aox > 2, i.e. it could be regarded as ‘X-ray 
weak’ compared to the average < aox >~ 16 for radio-quiet 
quasars (Yuan et al. 1998). In an ASCA observation the object 
was found in a high state, with a photon index of F = 1.98 and 
a Fe Ka line with EW ~ 300 eV. On shorter time scales, the 
X-ray flux in the 2-10 keV band is highly variable; the fastest 
variation detected is 60% in less than 2 x 10^^ s (Wang et al. 
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2000) . Historic light curves in the X-ray, UV, and optical bands 
indicate that the variability amplitude in the UV and optical 
(up to 70%) is much smaller than in X-rays (up to a factor 
of 10). The optical micro-variability of this object can actu¬ 
ally be driven by re-processing of the variable X-ray flux if 
only half of the absorbed X-rays are re-radiated in the optical- 
to-UV band. The comparison with Einstein (Kriss 1988) and 
EXOSAT (Malaguti et al. 1994) observations further showed 
that the quasar can be classihed as X-ray weak only in one out 
of five X-ray observations (for details see Wang et al. 2000). 

To what extent the intensity variations are related to spectral 
changes is difficult to examine in detail due to the poor statistics 
of the pointed ROSAT/ASCA spectra, especially if significant 
X-ray scattering is present. Corbin & Boroson (1996) mention 
that the quasar shows ‘associated absorption’ in an HST spec¬ 
trum and Wang et al. (2000) claim the existence of narrow as¬ 
sociated Lya absorption. There are no HST data for the C IV 
region, but co-added lUE data suggest C IV absorption consis¬ 
tent in velocity with the Lya absorption. 

2. The XMM-Newton observation 

PG 0844-1-349 was observed by XMM-Newton on November 
5, 2000. The EPIC PN was operated in Pull Window mode with 
medium biter with a total exposure of ^ 21 ksec. Both MOS 
cameras were in Large Window mode with medium biters as 
well, each yielding ^ 23.5 ksec of data. The two RGS chains 
were operated in Spectroscopy mode, both with exposures of 
~ 26 ksec. The OM camera was operated in Imaging mode 
and bve images were taken with the U-band biter (300-400 
nm). Por the data analysis we used the XMMSAS version 5.2.0 
for EPIC and OM, and version 5.3.0 for RGS. 

2.1. The light curve 

We determined the 0.5 — 10 keV PN light curve of 
PG 08444-349 by accumulating the photons in 400 sec bins 
from a circular region of radius r = 100" around the quasar po¬ 
sition; for this we selected single and double events with qual¬ 
ity bag = 0 (for details of the EPIC detectors see Ehle et al. 

2001) . Por the MOS data we took singles and doubles (Pattern 
0—12) as well and selected the photons from a circular region 
of radius r = 60" around the source. 

During the observation the background was rather noisy 
with some strong bares. The total background count rate gen¬ 
erally contributes only about 4% of the source count rate. 
However, in two bares the background count rate at energies 
above 3 keV reaches up to 40% of the source count rate and 
care has to be taken in the spectral analysis. In particular, we 
ignore data from the last 2 ksec of the observation as the back¬ 
ground count rate becomes very variable then. Since there are 
no differences between PN and MOS light curves, in order to 
increase the signal-to-noise ratio, we summed up the PN count 
rates (average 6.61 cts/s) and the two MOS count rates (av¬ 
erage MOS 1; 1.93 cts/s; MOS 2: 1.94 cts/s). The summed, 
background subtracted, EPIC light curve shown in Pig. starts 
at the switch on of the PN camera. There are about 2.5 ksec of 
MOS data previous to that, where the source stays on the same 



Fig. 1. Background subtracted summed PN + MOS light curve 
of PG 0844H-349; the time binning is 400 sec. 


high, bat count rate level. The RGS observation started around 
the same time as the PN, and lasted some 1.8 ksec longer, well 
into the period of baring background. 

The EPIC count rate (see Pig. starts to decrease after 
about 4 ksec, recovers for 2 ksec and then changes by about 
25% to a lower intensity level near the middle of the obser¬ 
vation. The bux drops occur in a smooth, nearly linear fash¬ 
ion at a rate of ~ 6.1 x 10“"* cts s“^. We will show in the 
next section that these changes are achromatic to high statisti¬ 
cal signibcance. With the spectral information obtained below, 
the K-corrected 0.5-10 keV luminosity during the brst part of 
the observation is ~ 2.5 x 10^"* erg s“^. We can use the ~ 
20% decline of the count rate over ~ 3500 s near the middle 
of the observation to estimate the lower limit of the radiative 
efficiency, p > 5 x AL / /S.t (Pabian 1984), and we ob¬ 

tain rj > 7 X 10“^ which is far below the theoretical limit for a 
Schwarzschild black hole. 

The OM provided U-blter images integrated for 4000 s 
each. The brightness of the source, obtained from the Pipeline 
Processing System (PPS), stayed approximately constant at 
14.00 mag; the largest deviation from this average value oc- 
cured during the period ^ 2-6 ksec after the start of the PN ex¬ 
posure, during the brst intensity dip, when the source was 0.01 
mag fainter than the average. This is, however, only a ^ 1.5a 
effect which implies that any variations of the optical bux are 
negligible compared to the changes in the X-ray band. 

2.2. Spectral properties 

Por the EPIC spectral analysis we used the latest available re¬ 
sponse matrices (version 6.1) issued in December 2001. Por the 
PN camera we exbacted single and double events with quality 
bag = 0 from a circular region with a radius of 45" around the 
source position. This radius includes ^ 88% of the source pho¬ 
tons (Ghizzardi & Molendi 2001) but avoids the gap between 
the detector chips. The background was taken on the same chip 
at a similar offset position. With a count rate of > 7 cts s“^ in 
the high state the PN detector, operated in Pull Window mode. 
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Fig. 2. Broken power law fit to the total PN + combined MOS 
data of PG 0844+349; the PN spectrum is the upper curve. 


showed strong indications of pile-up, clearly apparent from the 
XMMSAS task epatplot. We therefore discarded photons from 
the innermost 3x3 RAW pixels at the core of the point spread 
function from the spectral analysis. As a consequence, the nor¬ 
malizations of the models in the PN fits given in Table 1 are 
underestimating the actual photon flux. Further, we discarded 
the data of the last ~ 2 ksec of the observation, when the back¬ 
ground was strongly flaring, and we excluded times of high 
background, i.e. greater than 50 cts/s, from the Good Time 
Intervals (GTI) file. However, this last selection criterion influ¬ 
enced the quality of the fits only marginally. For the two MOS 
cameras we selected events with pattern < 12 and flag = 0 and 
a similar extraction radius as for the PN; the background was 
taken from an area on the same chips near the source. Pile-up 
did not influence the spectral fits of the MOS data noticeably; 
therefore the full data were used. 

The RGS spectra of PG 0844-1-349 were accumulated by 
selecting counts from a region centered on the source position 
and enclosing 90% of the point spread function in the cross 
dispersion direction, and making the first order selection by 
including 90% of the expected CCD pulse height distribution 
of the source photons. The background spectrum was obtained 
from two spatial regions excluding 95% of the source point 
spread function in the cross dispersion direction and perform¬ 
ing the same order selection as for the source. 

While the ASCA high state data could be well fitted by 
a simple absorbed power law (Wang et al. 2000), this model 
does not provide an acceptable fit to the PN data from the first, 
high state part of the observation. With a reduced = 1-42 
the fitted power law index is F = 2.57 and the residuals 
clearly show that at higher energies the spectrum gets flatter. 
A broken power law, fixing the absorption at the Galactic value 
Nh = 3.32 X 10^°cm“^ (Lockman & Savage 1995), gives an 
acceptable fit to the PN data with a reduced = 1-03 for 
372 d.o.f. (see Table 1). Fitting the same time interval with the 
MOS 1 and the MOS 2 data separately yielded similar accept¬ 
able fits. The fit parameters in Table 1 indicate that the two 
MOS spectra have nearly identical power law slopes. In both 
MOS cameras the residuals at the lowest energies show an S- 


shaped pattern with a maximum at ~ 0.5 keV, while the PN fits, 
with generally slightly flatter power laws and a higher break 
energy, show systematic positive residuals at ~ 0.6 keV. 

Applying the same model to the low state, i.e. the period be¬ 
tween 10 and 19 ksec after the start of the observation, gave an 
acceptable fit with nearly identical parameters, with the same 
‘typical’ differences between the three detectors as above (see 
Table 1). 

We repeated the fit over several other, shorter time intervals 
during the observation and always obtained acceptable fits with 
very similar parameters. We thus conclude that no statistically 
significant spectral variations occur during the flux changes of 
the source. Therefore for the subsequent fits we used the data 
from the full observation, again excluding the last 2 ksec and 
times of high background. 

Fig.|| shows the simultaneous broken power law fit to the 
whole data set for the PN and the combined MOS instruments. 
Due to the slight differences between the PN and MOS cam¬ 
eras, visible at ~ 0.5-0.6 keV and at higher energies, the qual¬ 
ity of the fit is only moderate (x^g^j = 1.11 for 1006 d.o.f.). 
Adding a Gaussian iron line does not improve the fit, mainly 
because of the low statistical significance of the data at higher 
energies. We therefore fitted the PN data separately for the 2-9 
keV range with a single power law and, secondly, with a power 
law plus a Gaussian line. We find a broad line at an energy 
of 6.25±0.13 keV (Ei^^est = 6.65 ± 0.13 keV in the quasar’s 
rest frame) and an equivalent width of ~ 235 eV. For a narrow 
line (i.e. fixing the Gaussian sigma at 10 eV) we get the same 
line energy and an equivalent width of ~ 120 eV. Including the 
lines improves the quality of the fits marginally (see Table 1), 
but only at a ~ 94% confidence level, according to an F-test. 
The small equivalent width found is atypical for NLSy 1 galax¬ 
ies, for which BeppoSax (Comastri 2000) and ASCA (Turner 
et al. 1998) generally find ionized iron lines with large equiva¬ 
lent widths (few 100 eV). 

The moderate quality of the combined PN and MOS fits 
and the fit residuals can be explained as resulting from the 
remaining uncertainties in the detector calibrations. However, 
they might also indicate that the chosen spectral model can be 
improved. 

In Fig. H we show the power law fit to the total PN data 
for the 2—10 keV range, extrapolated to lower energies. The 
slope of F ~ 2.04, consistent with ASCA, is probably not af¬ 
fected by contributions of more complex models for the low 
energy part of the spectrum. The ratio between the data and 
the model clearly shows the ‘gradual soft excess’ (Pounds & 
Reeves 2002) also found in other AGN with the high sensi¬ 
tivity and large bandwidth of XMM-Newton. While the hard 
power law tail seems to be well constrained, the exact shape 
of the spectrum at low energies remains uncertain. We there¬ 
fore tried several composite models, mainly to characterize the 
softer part of the spectrum. 

Models including a thermal component in addition to the 
hard power law fit the data quite well, as shown in Table 1; 
the relative contribution of the power law component can be 
inferred from the listed normalization. The power law + bbody 
fit provides two solutions: one in which the power law com¬ 
ponent represents the hard flux, the black body temperature is 
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Table 1. Results from spectral fitting assuming fixed Galactic Nh = 3.32 x 102Ocm-2. 


Period 

Detector 

Model 

r soft 

Ebreak 

(keV) 

r hard 

Power law norm. Xred/dof 

(10“® ph/keV/cm^/s) 

High 

PN 

bknpow 

2.71±0.02 

L63±0.12 

2.13±0.06 

L65±0.03 

1.03/372 


MOS 1 

bknpow 

2.81±0.04 

L02±0.08 

2.31±0.04 

2.96±0.12 

1.03/148 


MOS 2 

bknpow 

2.85±0.04 

L03±0.08 

2.31±0.04 

3.06±0.12 

1.10/145 

Low 

PN 

bknpow 

2.75±0.02 

L81±0.11 

2.08±0.06 

L32±0.02 

1.04/372 


MOS 1 

bknpow 

2.77±0.03 

L26±0.08 

2.27±0.04 

2.49±0.06 

0.92/167 


MOS 2 

bknpow 

2.93±0.03 

L05±0.06 

2.34±0.03 

2.45±0.08 

1.02/166 

Total 

PN 

bknpow 

2.72±0.01 

L80±0.08 

2.09±0.04 

L46±0.02 

0.98/537 

PN + MOS 1+2 bknpow 

2.77±0.09 

L33±0.03 

2.25±0.02 

(a) 

1.11/1006 

Total^i) 

PN 

pow 



2.05±0.05 

0.96±0.10 

0.83/188 


PN 

pow + gauss 

2.11±0.05 ( 

6.247±0.127(^> 

0.24±0.16(®^ 

0.96±0.10 

0.81/185 

Total 

PN 

pow + bbody 


L94±0.14H) 

2.73±0.01 

L45±0.02 

0.99/629 


PN 

pow + brems 


0.29±0.01H) 

2.11±0.03 

L02±0.06 

0.95/621 


PN 

pow + diskbb 


0.14±0.01H) 

2.23±0.02 

L17±0.04 

1.00/629 


PN 

pow + mekal 


0.26±0.01H) 

2.18±0.03 

L12±0.05 

0.95/627 


(a): Normalizations for the individual detectors allowed to be different. 

(1): Fit over the 2—9 keV energy band only; (2) Fitted line energy; (3): Line width (sigma) in keV; the photon flux in the line is 
(5.7±4.9) X 10“® photons/cm^/s. (4): Temperature of the additional model component in keV. 
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Fig. 3. Power law fit to the total PN data in the energy range 
2—10 keV; the fitted model is extrapolated to lower energies. 
The lower panel gives the ratio between data and model. 

similar to that of the diskbbody and another, perhaps less phys¬ 
ical solution, where the steep power law provides the soft flux, 
the hot bbody the hard flux (these parameters are given in Table 
1). All models are statistically acceptable (although the above 
mentioned systematic residuals at ~ 0.6 keV persist in all fits) 
and therefore the physical nature of the soft emission remains 
unclear. 

As the fitted slope of the spectrum at higher energies seems 
to harden with energy we tried to fit a curved continuum model 
(Fossati et al. 2000) but no good fit could be achieved. We fur¬ 
ther tried several composite cases fixing the hard power law 
index at the best fit slope of the 2—10 keV fit but we never 
obtained an acceptable fit: this seems to indicate that the hard 
power law is only the tail of a distribution and not a separate 
component. 
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Fig. 4. Power law + bremsstrahlung fit to the PN, RGS1 and 
RGS2 data of PG 0844+349. 

Nearly featureless spectra with a strong soft component are 
indicative of Comptonization of soft photons in the hot corona 
of an accretion disk (Haardt & Maraschi 1993). A fit with the 
Xspec Comptonization model comptt (Titarchuk 1994) failed 
to reproduce the hard power law and resulted in an unaccept¬ 
able fit (Xred = 1-97). Adding a hard power law to the comptt 
model (Xred = 1.09/1094 d.o.f) or fitting the sum of two 
Comptonization models with different temperatures and opti¬ 
cal depths of the scattering medium (Xred = 1.11/1093 d.o.f) 
resulted in fits nearly as good as the models in Table 1. For the 
double - comptt fit we either assumed the same temperature for 
the two soft photon components or left them free to vary inde¬ 
pendently. In all three cases the bulk of the flux up to ~ 3 keV 
originates from Comptonization of soft photons off a gas of kT 
~ 4.5 keV and optical depth r ~ 2.4. The high energy part of 
the spectrum is formed either by the power law component or 






W. Brinkmann et al.: XMM-Newton observations of PG 0844 + 349 


5 


the second Comptonization component with kT~ 16 keV and 
T ^ 3.3. The temperatures of the soft photons were htted to 
be kTo = 65^g5 eV in case of the extra power law compo¬ 
nent, kTo = 3l3^ eV (using the same values for the two soft 
components) and kTo(l) = 12132^ eV and kTo(2) = 2^2^^ 
eV for the case of two independent soft components, respec¬ 
tively. The temperatures of the Comptonizing electrons are at 
the low end of the expected temperature range for AGN (Haardt 
& Maraschi 1993) while the htted optical depths are higher 
than usually deduced. The temperatures of the soft photons ap¬ 
pear rather low; however, the parameters are very poorly con¬ 
strained because of the close coupling of temperature and opti¬ 
cal depth in the Comptonization models. Interestingly, the high 
optical depth implies that any rehection features (like the iron 
line) are suppressed by Compton scattering in the corona itself 
(Mattetal. 1997). 

The superior energy resolution of the RGS data could pro¬ 
vide an ultimate test for the nature of the soft emission. The 
same composite models tried on the PN data alone (see Table 
1) were htted simultaneously to the RGSl and RGS2 spec¬ 
tra, as well as to the PN. The hard power law slope and the 
temperature of the soft component were kept hx at the best ht 
value obtained from the hts of the PN alone, and the normal¬ 
izations for the three instruments were allowed to vary inde¬ 
pendently of each other (because of the PN data selection used 
to avoid pile-up). The low energy absorption column was hxed 
at the Galactic value. Values of x^^g/dof very similar to those 
in Table 1 for the PN alone were obtained (in the range 0.98- 
1.05/841), with the formally best ht model being the combina¬ 
tion of a power law and a bremsstrahlung component. This best 
ht and the data are shown in Fig.^ the ht of the pow/diskbb and 
pow/mekal are virtually indistinguishable from this. Although 
the signal-to-noise ratio of the RGS data is only moderate, it is 
clear that there is no evidence of emission or absorption struc¬ 
tures in the spectrum, indicating the absence of any additional 
source intrinsic absorbers. In particular, there is no evidence of 
excess emission in the RGS at ~0.6 keV which could explain 
the amplitude and the shape of the residuals observed in the PN 
ht. 

In the following, we use the broken power law model for 
the determination of the source’s energetics; however, it should 
be kept in mind that this is a rather artihcial representation of 
the source’s spectrum and that at the lowest energies the hux 
distribution might be quite different from a simple power law, 
introducing some uncertainty in the numerical values. 

The observed average un-absorbed hux in the high state in 
the 0.2-10.0 keV range amounts to 2.6xl0“^^erg cm“^s“^, 
while the average in the low state is 2.1xl0“^^erg cm“^s“^. 
With these Huxes the K-corrected X-ray luminosity of 

PG 0844-1-349 is in the range Lo, 2-10 keV = (4.1-4.9) x lO^^ 
erg s“^, using a Friedman cosmology with Hq = 50 km s“^ 
Mpc“^ and qo = 0.5. The average soft-band luminosity is 
Lo.5-2 keV ~ 6.4 X erg s“^, for the hard band we hnd 

L 2 -iokeV 6 X 10"^^ erg s“^. Thus PG 0844-1-349 is a 

quasar of rather low luminosity in the X-ray band. Using the 
monochromatic luminosity at 2 keV we derive an X-ray loud¬ 


ness Q!ox = —0.384 log(L 2 kev/T °) ~ 1-58, which is typ- 

ical for average radio-quiet quasars (Yuan et al. 1998). 

3. Discussion 

Using the normalizations of the power law hts the above 
analysis shows that PG 08444-349 was in a historically high 
state (see Table 1 of Wang et al. 2000) during the XMM- 
Newton observation, with an average 0.1 —2.4 keV luminosity 
of < 3 X 10^"* erg s Judging from the OM data the object 
was optically rather bright as well. An extrapolation of the soft 
power law into the optical band over-predicts the optical Hux 
by a factor of ~ 4, therefore the spectrum must break between 
the optical and the soft X-ray band. Nevertheless, most of the 
power is emitted in the soft X-ray / UV band. Taking the mass 
estimates for the central black hole from Kaspi et al. (2000) 
(2—3 X IO^Mq) this implies that PG 08444-349 must be a quite 
efficient accretor. The optical spectrum of PG 08444-349 re¬ 
sembles very much that of NLSy 1 galaxies (Boroson & Green 
1992 ) and thus an efficient accretion how as well as its X-ray 
variability would be in accordance with current models of these 
AGN. 

The X-ray spectrum of PG 08444-349 shows a relatively 
hat power law component at higher energies (Fhard ^ 2.25) 
and a strong, steep soft excess at energies < 2 keV. Both com¬ 
ponents are featureless and the inclusion of a broad or a narrow 
iron line around 6.7 keV improves the ht only marginally. The 
soft excess can be modeled with various components result¬ 
ing in hts with similar statistical significances. The RGS data 
do not indicate the presence of any spectral feature either. The 
excess of the soft Hux reaches a factor of ~2.5 over the extrap¬ 
olation of the hard power law. 

This kind of spectrum is not unusual amongst recent 
measurements of X-ray spectra of bright Seyferts / quasars 
with XMM-Newton (O’Brien et al. 2001, Page et al. 2001, 
Pounds & Reeves 2002). As a matter of fact, the spectrum 
of PG 08444-349 appears to be a carbon copy of that of the 
NLSyl galaxy PKS 0558-504 (Pounds & Reeves 2002, Fig. 
1), even with respect to the ‘big blue bump’ seen in both ob¬ 
jects (O'Brien et al. 2001). However, we also note that the 
power law slope of PG 08444-349 is slightly steeper than that 
of PKS 0558-504, which is radio-loud and shows correlations 
between X-ray brightness and hardness (Gliozzi et al. 2001), 
characteristics not shared by PG 08444-349 . 

The featureless spectra in these objects indicate that we are 
seeing the bare continuum disk emission from the quasars; thus 
disk Comptonization models, where the X-rays are produced 
via inverse Compton emission in a hot corona embedding a 
cooler accretion disk (e.g. Haardt & Maraschi 1993, Pounds 
et al. 1995), might provide a satisfactory physical explanation. 
Fits with available models for this scenario, however, yield pa¬ 
rameters for PG 08444-349 not typically found in other AGN, 
and not statistically preferred according to the above analysis; 
this might indicate that the actual physical conditions in the 
sources are more complex than our simplihed models can ac¬ 
count for. 

The large variety of possible Comptonization scenarios (for 
example, Haardt 1996, Zycki et al. 2001) does not allow a bet- 
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ter confinement of the parameters in the physical phase space. 
Stronger constraints might be deduceable from an analysis of 
the temporal behavior of the source. PG 0844+349 was known 
to be variable in X-rays from previous ROSAT and ASCA ob¬ 
servations, but only XMM-Newton is able to follow the flux 
changes on the shortest time scales. The X-ray flux (Fig. ba¬ 
sically changes during the observation from a higher to a lower 
level in a nearly linear fashion and the slopes of the intensity 
variations are very similar, indicating a rather well organized 
process. The light curve shows (for a NLSy 1) an atypically low 
variance and the intensity changes occur without any measur¬ 
able spectral changes in the X-ray band. 

The optical image of PG 0844-1-349 indicates that we are 
seeing the object nearly pole-on and therefore heavy obscura¬ 
tion of the X-ray emitting region appears unlikely. The emis¬ 
sion from a jet, changing its geometrical appearance, can very 
likely be ruled out as the source is dehnitely radio-quiet, un¬ 
less we are looking directly into a purely hydrodynamically 
outflowing jet for which there are no other observational in¬ 
dications. 

Achromatic flux changes can result from changes of the 
effective radiating area, ensuring that the spectral shape does 
not change noticeably. While the temperature of the cold mat¬ 
ter does not play a signihcant role in the models, changes 
in the optical depth r give rise to signihcant spectral vari¬ 
ability (Haardt et al. 1997). Most of the theoretical investi¬ 
gations are, however, concentrating on the higher energy part 
of the AGN spectrum (see, for example, Petrucci et al. 2001), 
and the parameter space of relevance here is only poorly ex¬ 
plored. In any case, changes of the radiating area are expected 
to happen on the dynamical time scale for Keplerian inhow, 
r ~ 9 10^ (Mbh/10^ Mq) sec, where is the 

Schwarzschild radius of the central black hole of mass Mbh- 
This estimate implies a rather low mass for the central black 
hole in PG 0844-1-349 and a compact emitting region. However, 
strong variability of Seyferts and Galactic black hole candi¬ 
dates indicates that the corona cannot be a uniform, continuous 
medium, unless it is geometrically thin (Celotti et al. 1992); 
moreover, observational evidence implies that the geometry of 
the coronal plasma cannot be slab like, but is made up of a 
number of distinct active regions (Haardt et al. 1994). 

In a popular class of models the corona is heated by mag¬ 
netic helds which rise up buoyantly from the disk where they 
reconnect and release their energy in flares. The energy stor¬ 
age in the corona is very likely the magnetic helds (Merloni & 
Fabian 2001) and the picture of the corona is that of a spread 
of active regions, of which only a few are large and dominate 
at any given time. The overall time scale of the evolution of the 
magnetic held conhguration is expected to be of the order of 
the Keplerian time scale again (Romanova et al. 1998) but in¬ 
dividual active regions can certainly evolve much more rapidly, 
depending of the magnetic held conhguration and the how con¬ 
ditions in the disk, which are largely unknown. 

4. Conclusions 

In a 20 ksec XMM-Newton observation the X-ray transient 
quasar PG 0844-1-349 was found in a historically high X-ray 


state with an average 0.2—lOkeV luminosity of Lq. 2-10 keV 
4.5 X 10''^'* erg s“^. During the observation the hux of 
PG 0844-1-349 varied achromatically by ^ 25% from a higher 
into a lower Hux state in a very smooth manner on a time scale 
of a few thousand seconds. 

The quasar showed a featureless spectrum which can be 
physically explained by Comptonization from hot electrons 
of the emission of an accretion disk. The rather low signal- 
to-noise ratio RGS spectra do not indicate the presence of 
any emission or absorption features either. The htted values 
of the slopes of the power law type spectrum at high ener¬ 
gies (Fhard 2.05) and those of the power law in the soft 
band (Fsoft ~ 2.75) differ slightly between the EPIC detec¬ 
tors and depend on the complexity of the models htted to 
the data. Current Comptonization models predict considerably 
lower temperatures and higher optical depths of the Compton 
scattering electrons than generally found in Seyfert galaxies 
and radio-quiet quasars (see e.g. Petrucci et al. 2001). As 
PG 0844-1-349 shares many of the characteristics of NLSyl 
galaxies and similar parameters were recently deduced for 
other NLSy 1 galaxies (O’ Brien et al. 2001, Page et al. 2001) we 
propose that Comptonization spectra with high optical depths 
and moderate plasma temperatures might be a distinguishing 
criterion for the accretion process in NLSyl galaxies. 

More sensitive XMM-Newton observations and a deeper 
theoretical exploration of this lower energy and lower temper¬ 
ature parameter space of Comptonization models might shed 
some new light into the physical properties of these objects. 
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Abstract. In a ~ 20 ksec XMM-Newton observation the X-ray transient radio-quiet quasar PG 0844-1-349 was found in 
a historically high state compared to previous X-ray observations. The quasar showed a featureless spectrum with a strong 
soft excess over the extrapolation of a hard power law. Comptonization models or a broken power law with Fsoft ~ 2.75, 
Fhard ~ 2.25 and a break energy of Ebreak ~ 1-35 keV represent acceptable descriptions of the spectral continuum. In the 
Comptonization models the temperature of the Comptonizing gas is considerably lower than generally found in (broad line) 
Seyfert galaxies whereas the optical depth is much higher. As a similar behavior has been seen in NLSyl galaxies, it might be 
an indicator of the different physical conditions in these two classes of AGN. During the XMM-Newton observation the flux of 
PG 0844-1-349 varied achromatically in a smooth, nearly linear fashion, hy ~ 25% on time scales of a few thousand seconds, 
which puts some constraints on current models of Comptonizing accretion disk coronae. 

Key words. Galaxies: active - quasars: individual: PG 0844-1-349 - X-rays: galaxies 


1. Introduction 

Quasars as a class of Active Galactic Nuclei (AGN) exhibit 
a vast diversity in their observed properties. Some of these, 
like the radio-loud - radio-quiet dichotomy, can be related to 
the physical characteristics of the accreting central black hole 
and to morphological differences of the host galaxies. Others 
can result from specific geometrical conditions of the accretion 
flow and from the particular viewing condition of the observer 
with respect to the quasar (see, for example, Elvis 2000). In 
many cases, however, the causes for the unusual properties of 
a quasar remain obscure and thus the study of objects display¬ 
ing atypical characteristics might yield important clues to the 
physical processes governing the quasar emission. 

One of the objects with several extraordinary properties 
is PG 08444-349, a well studied nearby (z=0.064), bright 
(mv=14.0), radio-quiet quasar from the Palomar Green sample 
(Schmidt & Green 1983). Deep radio observations reveal no 
nuclear source and the faint detection in the field (Si 49 ghz = 
0.3 mJy) might be related to an R = 20 mag object south-west 
of the quasar (Condon et al. 1987). PG 08444-349 is optically 
quite spectacular: it appears as a barred spiral with complex 
outer structure and a nearly equally bright companion galaxy 
(Hutchings & Crampton 1990). It belongs to the class of AGN 
with unusually strong Fe II emission (Wang et al. 1996) and ac¬ 
cording to the optical spectrum displayed in Boroson & Green 
(1992) the source shows all features of a Narrow-Line Seyfert 
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1 galaxy (NLSyl) - strong Fell emission and weak forbid¬ 
den narrow lines, although its FWHM(H/)) = 2420 km is 
slightly above the often used cut-off line of 2000 km (e.g. 
Osterbrock & Pogge, 1985). However, this cut-off for the clas¬ 
sification of NLSyl s is more artificial than physical, so we still 
might define this source as a NLSyl galaxy especially as it 
shows X-ray properties typical for this class of objects, i.e., a 
steep soft X-ray spectrum and strong variability (for reviews 
see Boiler et al. 2000). It is infrared-loud with a disturbed host 
galaxy (Clements 2000), and it is strongly variable at all wave¬ 
lengths on various time scales. The optical flux seems to have 
changed by ~70% between 1986 (Elvis et al. 1994) and 1993 
(Maoz et al. 1994) and intra-night variations with amplitudes of 
AB = 0.08 mag were detected by Jang et al. (1997). From re¬ 
verberation measurements Kaspi et al. (2000) deduced a mass 
of ~ 2.5 X 10^ Mq for the quasar. 

Compared to the ROSAT All Sky Survey (Yuan et al. 1998) 
the soft X-ray flux of PG 08444-349 had decreased by a factor 
of ~6 in a pointed ROSAT observation 6 months later, with¬ 
out any noticeable changes in the spectral slope (Rachen et al. 
1996). During the X-ray low state PG 08444-349 had an X- 
ray loudness of aox > 2, i.e. it could be regarded as ‘X-ray 
weak’ compared to the average < aox >~ 1.6 for radio-quiet 
quasars (Yuan et al. 1998). In an ASCA observation the object 
was found in a high state, with a photon index of L = 1.98 and 
a Fe Ka line with EW ~ 300 eV. On shorter time scales, the 
X-ray flux in the 2-10 keV band is highly variable; the fastest 
variation detected is 60% in less than 2 x 10^^ s (Wang et al. 
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2000) . Historic light curves in the X-ray, UV, and optical bands 
indicate that the variability amplitude in the UV and optical 
(up to 70%) is much smaller than in X-rays (up to a factor 
of 10). The optical micro-variability of this object can actu¬ 
ally be driven by re-processing of the variable X-ray flux if 
only half of the absorbed X-rays are re-radiated in the optical- 
to-UV band. The comparison with Einstein (Kriss 1988) and 
EXOSAT (Malaguti et al. 1994) observations further showed 
that the quasar can be classified as X-ray weak only in one out 
of five X-ray observations (for details see Wang et al. 2000). 

To what extent the intensity variations are related to spectral 
changes is difficult to examine in detail due to the poor statistics 
of the pointed ROSAT/ASCA spectra, especially if significant 
X-ray scattering is present. Corbin & Boroson (1996) mention 
that the quasar shows ‘associated absorption’ in an HST spec¬ 
trum and Wang et al. (2000) claim the existence of narrow as¬ 
sociated Lya absorption. There are no HST data for the C IV 
region, but co-added lUE data suggest C IV absorption consis¬ 
tent in velocity with the Lya absorption. 

2. The XMM-Newton observation 

PG 0844-1-349 was observed by XMM-Newton on November 
5,2000. The EPIC PN was operated in Pull Window mode with 
medium filter with a total exposure of ~ 21 ksec. Both MOS 
cameras were in Large Window mode with medium filters as 
well, each yielding ~ 23.5 ksec of data. The two RGS chains 
were operated in Spectroscopy mode, both with exposures of 
~ 26 ksec. The OM camera was operated in Imaging mode 
and five images were taken with the U-band filter (300^00 
nm). Por the data analysis we used the XMMSAS version 5.2.0 
for EPIC and OM, and version 5.3.0 for RGS. 

2.1. The light curve 

We determined the 0.5 — 10 keV PN light curve of 
PG 0844-1-349 by accumulating the photons in 400 sec bins 
from a circular region of radius r = 100 " around the quasar po¬ 
sition; for this we selected single and double events with qual¬ 
ity flag = 0 (for details of the EPIC detectors see Ehle et al. 

2001) . Por the MOS data we took singles and doubles (Pattern 
0 — 12 ) as well and selected the photons from a circular region 
of radius r = 60" around the source. 

During the observation the background was rather noisy 
with some strong flares. The total background count rate gen¬ 
erally contributes only about 4% of the source count rate. 
However, in two flares the background count rate at energies 
above 3 keV reaches up to 40% of the source count rate and 
care has to be taken in the spectral analysis. In particular, we 
ignore data from the last 2 ksec of the observation as the back¬ 
ground count rate becomes very variable then. Since there are 
no differences between PN and MOS light curves, in order to 
increase the signal-to-noise ratio, we summed up the PN count 
rates (average 6.61 cts/s) and the two MOS count rates (av¬ 
erage MOS 1: 1.93 cts/s; MOS 2: 1.94 cts/s). The summed, 
background subtracted, EPIC light curve shown in Pig. 1 starts 
at the switch on of the PN camera. There are about 2.5 ksec of 
MOS data previous to that, where the source stays on the same 
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Fig. 1. Background subtracted summed PN -l- MOS light curve of PG 
0844-1-349; the time binning is 400 sec. 

high, flat count rate level. The RGS observation started around 
the same time as the PN, and lasted some 1.8 ksec longer, well 
into the period of flaring background. 

The EPIC count rate (see Pig. 1) starts to decrease after 
about 4 ksec, recovers for ~ 2 ksec and then changes by about 
25% to a lower intensity level near the middle of the obser¬ 
vation. The flux drops occur in a smooth, nearly linear fash¬ 
ion at a rate of ~ 6.1 x 10“^^ cts s“^. We will show in the 
next section that these changes are achromatic to high statisti¬ 
cal significance. With the spectral information obtained below, 
the K-corrected 0.5-10 keV luminosity during the first part of 
the observation is ~ 2.5 x 10^^^ erg s^^. We can use the ~ 
20% decline of the count rate over ~ 3500 s near the middle 
of the observation to estimate the lower limit of the radiative 
efficiency, r] > 5 x lO^^^^AL/Af (Pabian 1984), and we ob¬ 
tain 7] > 7 X 10^^ which is far below the theoretical limit for a 
Schwarzschild black hole. 

The OM provided U-filter images integrated for ~ 4000 s 
each. The brightness of the source, obtained from the Pipeline 
Processing System (PPS), stayed approximately constant at 
14.00 mag; the largest deviation from this average value oc- 
cured during the period ~ 2-6 ksec after the start of the PN ex¬ 
posure, during the first intensity dip, when the source was 0.01 
mag fainter than the average. This is, however, only a ~ 1.5a 
effect which implies that any variations of the optical flux are 
negligible compared to the changes in the X-ray band. 

2.2. Spectral properties 

Por the EPIC spectral analysis we used the latest available re¬ 
sponse matrices (version 6.1) issued in December 2001. Por the 
PN camera we extracted single and double events with quality 
flag = 0 from a circular region with a radius of 45" around the 
source position. This radius includes ~ 88% of the source pho¬ 
tons (Ghizzardi & Molendi 2001) but avoids the gap between 
the detector chips. The background was taken on the same chip 
at a similar offset position. With a count rate of > 7 cts s^^ in 
the high state the PN detector, operated in Pull Window mode, 
showed strong indications of pile-up, clearly apparent from the 
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Fig. 2. Broken power law fit to the total PN + combined MOS data of 
PG 0844+349; the PN spectrum is the upper curve. 


XMMSAS task epatplot. We therefore discarded photons from 
the innermost 3x3 RAW pixels at the core of the point spread 
function from the spectral analysis. As a consequence, the nor¬ 
malizations of the models in the PN fits given in Table 1 are 
underestimating the actual photon flux. Further, we discarded 
the data of the last ~ 2 ksec of the observation, when the back¬ 
ground was strongly flaring, and we excluded times of high 
background, i.e. greater than 50 cts/s, from the Good Time 
Intervals (GTI) file. However, this last selection criterion influ¬ 
enced the quality of the fits only marginally. For the two MOS 
cameras we selected events with pattern <12 and flag = 0 and 
a similar extraction radius as for the PN; the background was 
taken from an area on the same chips near the source. Pile-up 
did not influence the spectral fits of the MOS data noticeably; 
therefore the full data were used. 

The RGS spectra of PG 0844-1-349 were accumulated by 
selecting counts from a region centered on the source position 
and enclosing 90% of the point spread function in the cross 
dispersion direction, and making the first order selection by 
including 90% of the expected CCD pulse height distribution 
of the source photons. The background spectrum was obtained 
from two spatial regions excluding 95% of the source point 
spread function in the cross dispersion direction and perform¬ 
ing the same order selection as for the source. 

While the ASCA high state data could be well fitted by 
a simple absorbed power law (Wang et al. 2000), this model 
does not provide an acceptable fit to the PN data from the first, 
high state part of the observation. With a reduced = 1-42 
the fitted power law index is F = 2.57 and the residuals 
clearly show that at higher energies the spectrum gets flatter. 
A broken power law, fixing the absorption at the Galactic value 
Nh = 3.32 X 10^°cm“^ (Lockman & Savage 1995), gives an 
acceptable fit to the PN data with a reduced = 103 for 
372 d.o.f. (see Table 1). Fitting the same time interval with the 
MOS 1 and the MOS 2 data separately yielded similar accept¬ 
able fits. The fit parameters in Table 1 indicate that the two 
MOS spectra have nearly identical power law slopes. In both 
MOS cameras the residuals at the lowest energies show an S- 
shaped pattern with a maximum at ~ 0.5 keV, while the PN fits. 


with generally slightly flatter power laws and a higher break 
energy, show systematic positive residuals at ~ 0.6 keV. 

Applying the same model to the low state, i.e. the period be¬ 
tween 10 and 19 ksec after the start of the observation, gave an 
acceptable fit with nearly identical parameters, with the same 
‘typical’ differences between the three detectors as above (see 
Table 1). 

We repeated the fit over several other, shorter time intervals 
during the observation and always obtained acceptable fits with 
very similar parameters. We thus conclude that no statistically 
significant spectral variations occur during the flux changes of 
the source. Therefore for the subsequent fits we used the data 
from the full observation, again excluding the last 2 ksec and 
times of high background. 

Fig.2 shows the simultaneous broken power law fit to the 
whole data set for the PN and the combined MOS instruments. 
Due to the slight differences between the PN and MOS cam¬ 
eras, visible at ~ 0.5-0.6 keV and at higher energies, the qual¬ 
ity of the fit is only moderate (Xred “ 1006 d.o.f.). 

Adding a Gaussian iron line does not improve the fit, mainly 
because of the low statistical significance of the data at higher 
energies. We therefore fitted the PN data separately for the 2-9 
keV range with a single power law and, secondly, with a power 
law plus a Gaussian line. We find a broad line at an energy 
of 6.25±0.13 keV (Ei rest = 6.65 ± 0.13 keV in the quasar’s 
rest frame) and an equivalent width of ~ 235 eV. For a narrow 
line (i.e. fixing the Gaussian sigma at 10 eV) we get the same 
line energy and an equivalent width of ~ 120 eV. Including the 
lines improves the quality of the fits marginally (see Table 1), 
but only at a ~ 94% confidence level, according to an F-test. 
The small equivalent width found is atypical for NLSyl galax¬ 
ies, for which BeppoSax (Comastri 2000) and ASCA (Turner 
et al. 1998) generally find ionized iron lines with large equiva¬ 
lent widths (few 100 eV). 

The moderate quality of the combined PN and MOS fits 
and the fit residuals can be explained as resulting from the 
remaining uncertainties in the detector calibrations. However, 
they might also indicate that the chosen spectral model can be 
improved. 

In Fig. 3 we show the power law fit to the total PN data 
for the 2—10 keV range, extrapolated to lower energies. The 
slope of F ~ 2.04, consistent with ASCA, is probably not af¬ 
fected by contributions of more complex models for the low 
energy part of the spectrum. The ratio between the data and 
the model clearly shows the ‘gradual soft excess’ (Pounds & 
Reeves 2002) also found in other AGN with the high sensi¬ 
tivity and large bandwidth of XMM-Newton. While the hard 
power law tail seems to be well constrained, the exact shape 
of the spectrum at low energies remains uncertain. We there¬ 
fore tried several composite models, mainly to characterize the 
softer part of the spectrum. 

Models including a thermal component in addition to the 
hard power law fit the data quite well, as shown in Table 1; 
the relative contribution of the power law component can be 
inferred from the listed normalization. The power law + bbody 
fit provides two solutions: one in which the power law com¬ 
ponent represents the hard flux, the black body temperature is 
similar to that of the diskbbody and another, perhaps less phys- 
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Table 1. Results from spectral fitting assuming fixed Galactic Nh = 3.32 x 10^°cm 


Period 

Detector 

Model 

r soft 

Ebreak 

(keV) 

E hard 

Power law norm. xJed^dof 

(10^® ph/keV/cm^/s) 

High 

PN 

hknpow 

2.71±0.02 

L63±0.12 

2.13±0.06 

L65±0.03 

1.03/372 


MOS 1 

bknpow 

2.81±0.04 

L02±0.08 

2.31±0.04 

2.96±0.12 

1.03/148 


MOS 2 

bknpow 

2.85±0.04 

L03±0.08 

2.31±0.04 

3.06±0.12 

1.10/145 

Low 

PN 

bknpow 

2.75±0.02 

L81±0.11 

2.08±0.06 

L32±0.02 

1.04/372 


MOS 1 

bknpow 

2.77±0.03 

L26±0.08 

2.27±0.04 

2.49±0.06 

0.92/167 


MOS 2 

bknpow 

2.93±0.03 

L05±0.06 

2.34±0.03 

2.45±0.08 

1.02/166 

Total 

PN 

bknpow 

2.72±0.01 

L80±0.08 

2.09±0.04 

L46±0.02 

0.98/537 

PN + MOS 1+2 bknpow 

2.77±0.09 

L33±0.03 

2.25±0.02 

(a) 

1.11/1006 

Totak^) 

PN 

pow 



2.05±0.05 

0.96±0.10 

0.83/188 


PN 

pow + gauss 

2.11±0.05 

6.247±0.127(^> 

0.24±0.16(®> 

0.96±0.10 

0.81/185 

Total 

PN 

pow + bbody 


L94±0.14('*) 

2.73±0.01 

L45±0.02 

0.99/629 


PN 

pow + brems 


0.29±0.0l('*) 

2.11±0.03 

L02±0.06 

0.95/621 


PN 

pow + diskbb 


0.14±0.01H) 

2.23±0.02 

L17±0.04 

1.00/629 


PN 

pow + mekal 


0.26±0.01H) 

2.18±0.03 

L12±0.05 

0.95/627 


(a): Normalizations for the individual detectors allowed to be different. 

(1): Fit over the 2—9 keV energy band only; (2) Fitted line energy; (3): Line width (sigma) in keV; the photon flux in the line is 
(5.7±4.9) X 10^® photons/cm^/s. (4): Temperature of the additional model component in keV. 
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Fig. 3. Power law fit to the total PN data in the energy range 2—10 keV; 
the fitted model is extrapolated to lower energies. The lower panel 
gives the ratio between data and model. 


ical solution, where the steep power law provides the soft flux, 
the hot bbody the hard flux (these parameters are given in Table 
1). All models are statistically acceptable (although the above 
mentioned systematic residuals at ~ 0.6 keV persist in all fits) 
and therefore the physical nature of the soft emission remains 
unclear. 

As the fitted slope of the spectrum at higher energies seems 
to harden with energy we tried to fit a curved continuum model 
(Fossati et al. 2000) but no good fit could be achieved. We fur¬ 
ther tried several composite cases fixing the hard power law 
index at the best fit slope of the 2—10 keV fit but we never 
obtained an acceptable fit: this seems to indicate that the hard 
power law is only the tail of a distribution and not a separate 
component. 
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Fig. 4. Power law + hremsstrahlung fit to the PN, RGSl and RGS2 
data of PG 0844+349. 

Nearly featureless spectra with a strong soft component are 
indicative of Comptonization of soft photons in the hot corona 
of an accretion disk (Haardt & Maraschi 1993). A fit with the 
Xspec Comptonization model comptt (Titarchuk 1994) failed 
to reproduce the hard power law and resulted in an unaccept¬ 
able fit iXied — 1-97). Adding a hard power law to the comptt 
model (x^gjj = 1.09/1094 d.o.f) or fitting the sum of two 
Comptonization models with different temperatures and opti¬ 
cal depths of the scattering medium (Xied ~ 1.11/1093 d.o.f) 
resulted in fits nearly as good as the models in Table 1. For the 
double - comptt fit we either assumed the same temperature for 
the two soft photon components or left them free to vary inde¬ 
pendently. In all three cases the bulk of the flux up to ~ 3 keV 
originates from Comptonization of soft photons off a gas of kT 
~ 4.5 keV and optical depth r ~ 2.4. The high energy part of 
the spectrum is formed either by the power law component or 
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the second Comptonization component with kT~ 16 keV and 
T ~ 3.3. The temperatures of the soft photons were fitted to 
be kTo = 65^05 eV in case of the extra power law compo¬ 
nent, kTo = 3l3^ eV (using the same values for the two soft 
components) and kTo(l) = 121^®^ eV and kTo(2) = 212^^^ 
eV for the case of two independent soft components, respec¬ 
tively. The temperatures of the Comptonizing electrons are at 
the low end of the expected temperature range for AGN (Haardt 
& Maraschi 1993) while the fitted optical depths are higher 
than usually deduced. The temperatures of the soft photons ap¬ 
pear rather low; however, the parameters are very poorly con¬ 
strained because of the close coupling of temperature and opti¬ 
cal depth in the Comptonization models. Interestingly, the high 
optical depth implies that any reflection features (like the iron 
line) are suppressed by Compton scattering in the corona itself 
(Mattetal. 1997). 

The superior energy resolution of the RGS data could pro¬ 
vide an ultimate test for the nature of the soft emission. The 
same composite models tried on the PN data alone (see Table 
1) were fitted simultaneously to the RGSl and RGS2 spec¬ 
tra, as well as to the PN. The hard power law slope and the 
temperature of the soft component were kept fix at the best fit 
value obtained from the fits of the PN alone, and the normal¬ 
izations for the three instruments were allowed to vary inde¬ 
pendently of each other (because of the PN data selection used 
to avoid pile-up). The low energy absorption column was fixed 
at the Galactic value. Values of x^^^/dof very similar to those 
in Table 1 for the PN alone were obtained (in the range 0.98- 
1.05/841), with the formally best fit model being the combina¬ 
tion of a power law and a bremsstrahlung component. This best 
fit and the data are shown in Fig.4; the fit of the pow/diskbb and 
pow/mekal are virtually indistinguishable from this. Although 
the signal-to-noise ratio of the RGS data is only moderate, it is 
clear that there is no evidence of emission or absorption struc¬ 
tures in the spectrum, indicating the absence of any additional 
source intrinsic absorbers. In particular, there is no evidence of 
excess emission in the RGS at ~0.6 keV which could explain 
the amplitude and the shape of the residuals observed in the PN 
fit. 

In the following, we use the broken power law model for 
the determination of the source’s energetics; however, it should 
be kept in mind that this is a rather artificial representation of 
the source’s spectrum and that at the lowest energies the flux 
distribution might be quite different from a simple power law, 
introducing some uncertainty in the numerical values. 

The observed average un-absorbed flux in the high state in 
the 0.2-10.0 keV range amounts to 2.6xl0^^^erg cm^^s^^, 
while the average in the low state is 2.1xl0^^^erg cm^^s^^. 
With these fluxes the K-corrected X-ray luminosity of 
PG 08444-349 is in the range Lo. 2-10 keV = (4.1 - 4.9) x 10^^ 
erg s“^, using a Friedman cosmology with Hq = 50 km s“^ 
Mpc^^ and qo = 0.5. The average soft-band luminosity is 
Lo. 5-2 keV ~ 6.4 X 10"^^ erg s^^, for the hard band we find 
L 2 -i 0 keV ~ 6 X 10^3 erg s-T Thus PG 08444-349 is a 
quasar of rather low luminosity in the X-ray band. Using the 
monochromatic luminosity at 2 keV we derive an X-ray loud¬ 


ness Q!ox = —0.384 log(L 2 kev/T °) ~ 1-58, which is typ- 

2 500 ^4 

ical for average radio-quiet quasars (Yuan et al. 1998). 

3. Discussion 

Using the normalizations of the power law fits the above 
analysis shows that PG 08444-349 was in a historically high 
state (see Table 1 of Wang et al. 2000) during the XMM- 
Newton observation, with an average 0.1 —2.4 keV luminosity 
of < 3 X 10"*"^ erg s“^. Judging from the OM data the object 
was optically rather bright as well. An extrapolation of the soft 
power law into the optical band over-predicts the optical flux 
by a factor of ~ 4, therefore the spectrum must break between 
the optical and the soft X-ray band. Nevertheless, most of the 
power is emitted in the soft X-ray / UV band. Taking the mass 
estimates for the central black hole from Kaspi et al. (2000) 
(2—3 X IO^Mq) this implies that PG 08444-349 must be a quite 
efficient accretor. The optical spectrum of PG 08444-349 re¬ 
sembles very much that of NLSyl galaxies (Boroson & Green 
1992 ) and thus an efficient accretion flow as well as its X-ray 
variability would be in accordance with current models of these 
AGN. 

The X-ray spectrum of PG 08444-349 shows a relatively 
flat power law component at higher energies (Fhard ~ 2.25) 
and a strong, steep soft excess at energies < 2 keV. Both com¬ 
ponents are featureless and the inclusion of a broad or a narrow 
iron line around 6.7 keV improves the fit only marginally. The 
soft excess can be modeled with various components result¬ 
ing in fits with similar statistical significances. The RGS data 
do not indicate the presence of any spectral feature either. The 
excess of the soft flux reaches a factor of ~2.5 over the extrap¬ 
olation of the hard power law. 

This kind of spectrum is not unusual amongst recent 
measurements of X-ray spectra of bright Seyferts / quasars 
with XMM-Newton (O’Brien et al. 2001, Page et al. 2001, 
Pounds & Reeves 2002). As a matter of fact, the spectrum 
of PG 08444-349 appears to be a carbon copy of that of the 
NLSyl galaxy PKS 0558-504 (Pounds & Reeves 2002, Fig. 
1 ), even with respect to the ‘big blue bump’ seen in both ob¬ 
jects (O'Brien et al. 2001). However, we also note that the 
power law slope of PG 08444-349 is slightly steeper than that 
of PKS 0558-504, which is radio-loud and shows correlations 
between X-ray brightness and hardness (Gliozzi et al. 2001), 
characteristics not shared by PG 08444-349 . 

The featureless spectra in these objects indicate that we are 
seeing the bare continuum disk emission from the quasars: thus 
disk Comptonization models, where the X-rays are produced 
via inverse Compton emission in a hot corona embedding a 
cooler accretion disk (e.g. Haardt & Maraschi 1993, Pounds 
et al. 1995), might provide a satisfactory physical explanation. 
Fits with available models for this scenario, however, yield pa¬ 
rameters for PG 08444-349 not typically found in other AGN, 
and not statistically preferred according to the above analysis: 
this might indicate that the actual physical conditions in the 
sources are more complex than our simplified models can ac¬ 
count for. 

The large variety of possible Comptonization scenarios (for 
example, Haardt 1996, Zycki et al. 2001) does not allow a bet- 
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ter confinement of the parameters in the physical phase space. 
Stronger constraints might be deduceable from an analysis of 
the temporal behavior of the source. PG 0844+349 was known 
to be variable in X-rays from previous ROSAT and ASCA ob¬ 
servations, but only XMM-Newton is able to follow the flux 
changes on the shortest time scales. The X-ray flux (Fig. 1) ba¬ 
sically changes during the observation from a higher to a lower 
level in a nearly linear fashion and the slopes of the intensity 
variations are very similar, indicating a rather well organized 
process. The light curve shows (for a NLSy 1) an atypically low 
variance and the intensity changes occur without any measur¬ 
able spectral changes in the X-ray band. 

The optical image of PG 0844-1-349 indicates that we are 
seeing the object nearly pole-on and therefore heavy obscura¬ 
tion of the X-ray emitting region appears unlikely. The emis¬ 
sion from a jet, changing its geometrical appearance, can very 
likely be ruled out as the source is definitely radio-quiet, un¬ 
less we are looking directly into a purely hydrodynamically 
outflowing jet for which there are no other observational in¬ 
dications. 

Achromatic flux changes can result from changes of the 
effective radiating area, ensuring that the spectral shape does 
not change noticeably. While the temperature of the cold mat¬ 
ter does not play a significant role in the models, changes 
in the optical depth r give rise to significant spectral vari¬ 
ability (Haardt et al. 1997). Most of the theoretical investi¬ 
gations are, however, concentrating on the higher energy part 
of the AGN spectrum (see, for example, Petrucci et al. 2001), 
and the parameter space of relevance here is only poorly ex¬ 
plored. In any case, changes of the radiating area are expected 
to happen on the dynamical time scale for Keplerian inflow, 
r ~ 9 10^ (r/Rg)^/^ (Mbh/10^ 1^©) where Rg is the 
Schwarzschild radius of the central black hole of mass Mbh- 
This estimate implies a rather low mass for the central black 
hole in PG 0844-1-349 and a compact emitting region. However, 
strong variability of Seyferts and Galactic black hole candi¬ 
dates indicates that the corona cannot be a uniform, continuous 
medium, unless it is geometrically thin (Celotti et al. 1992); 
moreover, observational evidence implies that the geometry of 
the coronal plasma cannot be slab like, but is made up of a 
number of distinct active regions (Haardt et al. 1994). 

In a popular class of models the corona is heated by mag¬ 
netic fields which rise up buoyantly from the disk where they 
reconnect and release their energy in flares. The energy stor¬ 
age in the corona is very likely the magnetic fields (Merloni & 
Fabian 2001) and the picture of the corona is that of a spread 
of active regions, of which only a few are large and dominate 
at any given time. The overall time scale of the evolution of the 
magnetic field configuration is expected to be of the order of 
the Keplerian time scale again (Romanova et al. 1998) but in¬ 
dividual active regions can certainly evolve much more rapidly, 
depending of the magnetic field configuration and the flow con¬ 
ditions in the disk, which are largely unknown. 

4. Conclusions 

In a 20 ksec XMM-Newton observation the X-ray transient 
quasar PG 0844-1-349 was found in a historically high X-ray 


state with an average 0.2—10 keV luminosity of Lo, 2 -io kev ~ 
4.5 X 10"^^ erg s^^. During the observation the flux of 
PG 0844-1-349 varied achromatically by ~ 25% from a higher 
into a lower flux state in a very smooth manner on a time scale 
of a few thousand seconds. 

The quasar showed a featureless spectrum which can be 
physically explained by Comptonization from hot electrons 
of the emission of an accretion disk. The rather low signal- 
to-noise ratio RGS spectra do not indicate the presence of 
any emission or absorption features either. The fitted values 
of the slopes of the power law type spectrum at high ener¬ 
gies (Fhard ~ 2.05) and those of the power law in the soft 
band (Fgoft ~ 2.75) differ slightly between the EPIC detec¬ 
tors and depend on the complexity of the models fitted to 
the data. Current Comptonization models predict considerably 
lower temperatures and higher optical depths of the Compton 
scattering electrons than generally found in Seyfert galaxies 
and radio-quiet quasars (see e.g. Petrucci et al. 2001). As 
PG 0844-1-349 shares many of the characteristics of NLSyl 
galaxies and similar parameters were recently deduced for 
other NLSyl galaxies (O’Brien et al. 2001, Page et al. 2001) we 
propose that Comptonization spectra with high optical depths 
and moderate plasma temperatures might be a distinguishing 
criterion for the accretion process in NLSyl galaxies. 

More sensitive XMM-Newton observations and a deeper 
theoretical exploration of this lower energy and lower temper¬ 
ature parameter space of Comptonization models might shed 
some new light into the physical properties of these objects. 
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